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We propose that the leptonic cosmic ray signals seen by PAMELA and ATIC result from 
the annihilation or decay of dark matter particles via states of a leptonic Higgs doublet to r 
leptons, linking cosmic ray signals of dark matter to LHC signals of the Higgs sector. The 
states of the leptonic Higgs doublet are lighter than about 200 GeV, yielding large fr and 
TTTT event rates at the LHC. Simple models are given for the dark matter particle and its 
interactions with the leptonic Higgs, for cosmic ray signals arising from both annihilations 
and decays in the galactic halo. For the case of annihilations, cosmic photon and neutrino 
signals are on the verge of discovery. 
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I. INTRODUCTION 

Recent observations of high-energy electron and positron cosmic ray spectra have generated tremen- 
dous interest, as they might provide the first non-gravitational evidence for Dark Matter (DM). The 
PAMELA [l| experiment reports an excess of positrons in the few GeV to 100 GeV range, providing 
further support to the earlier results of HEAT [2] and AMS 3]. In addition, results from the ATIC [J] 
and PPB-BETS 0] balloon experiments suggest an excess of electrons and positrons in the 300 GeV to 
600 GeV range. 

While these observations have conventional astrophysical interpretations, they may result from anni- 
hilations or decays of DM particles in the galactic halo. Indeed, the PAMELA and ATIC data reinforce 
each other, since, for a certain range of DM masses, they have a unified interpretation. However, DM 
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explanations for the leptonic cosmic ray excesses face two interesting challenges. First, for annihilating 
DM these signals require that the annihilation cross-section for DM particles is typically two to three 
orders of magnitude larger than that expected from the thermal freezeout of WIMP DM. On the other 
hand, for decaying DM, the life-time of the DM particles must be extremely large, of O(10^^~^^) sec- 
onds. Second, the signals apparently require annihilations or decays dominantly into leptons rather than 
hadrons, since there is no reported excess in anti-proton cosmic rays. Inspite of these challenges, many 
sapers with different models of DM have already appeared in the literature, utilizing both annihilations 



papers witn diirerent models oi JJJ 
H 0, B B and decays Q, Q, Q • 



For annihilating DM, many of these models try to explain the required large annihilation cross-section 
by a Sommerfeld enhancement which is operative at the non-relativistic velocities (/3 ~ 10 ) in 
the galactic halo, while still having a standard thermal relic abundance applicable at the time of DM 
freeze-out. There are several possibilities for understanding why the products of the annihilation are 
dominantly leptonic rather than hadronic. One possibility is kinematics: the annihilation products of the 
DM particles are not heavy enough to decay into quarks, gauge bosons, and Higgs bosons, which have a 
large hadronic branching ratio, and hence decay only to electrons and muons (and possibly taus). 

Another possibility is a symmetry, rather than kinematics, to understand why the DM annihilation 
or decay products are dominantly leptonic. In this work, we study a DM sector coupling to the visible 
sector through Higgs messengers which couple to leptons due to a symmetry. The plan of the paper is as 
follows. In the next section, we motivate this possibility from a general perspective, stressing that this is 
a natural implementation of the hypothesis that the DM is a WIMP, with mass and interactions broadly 
governed by the mass scale of weak interactions. In section IIIIl we give the yukawa interactions for the 
minimal leptonic Higgs theory, and write them in a mass eigenstate basis for both the Higgs and the 
quarks and leptons. There are two Higgs scalars, h and H, one pseudoscalar. A, and one charged Higgs 
boson H^, each with matter interactions that are determined by the ratio of vevs tan/3 and the Higgs 
mixing angle a. These interactions are quite unlike those of the usually considered two Higgs doublet 
model or of the MSSM. The LEP mass limits for h, H, A and H'^ are given, as well as constraints that 
follow from the cosmic ray signals. 

In section IIVI the possible classes of DM annihilations and decays through the leptonic Higgs states 
H, A and are discussed. The common feature is multi-r final states. For a particular annihilation 
channel the cosmic-ray electron and positrons signals are derived and compared to the PAMELA and 
ATIC data. In general, if the PAMELA and ATIC data result from enhanced galactic DM annihilation, 
then significant fluxes of photons Q Q Q and neutrinos 0. Q »e also expected^ On the one hand 
this could allow for a crucial confirmation of the DM nature of the signal, while on the other hand there 
is frequently some tension with present limits on photon and neutrino fluxes. We calculate photon fluxes 
in the leptonic Higgs model from DM annihilations in both the galactic center of the Milky Way and in 
the dwarf galaxy Sagittarius, and also discuss the neutrino flux. 

In section |Vl we study the implications for Higgs signals at the LHC. The LHC Higgs signals are 
necessarily correlated with the cosmic-ray signals, with the leptonic Higgs states H, A and decaying 
dominantly to final states involving the r lepton. The collider phenomenology is very rich, and quite 
unlike that of conventional two Higgs doublet models, such as the MSSM. In section IVTl we study some 
simple models for the DM particle and its couplings to the Higgs sector, including the cases that the 
DM particle is derived from electroweak singlet and doublet fermions or scalars. Since our dark matter 
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particle is heavy, in the few TeV region, for DM annihilations the coupling strength to the leptonic 
Higgs is typically quite strong. In addition, a Sommerfeld enhancement of the annihilation cross section 
typically results from the exchange of the leptonic Higgs, which is lighter than 2M\y. In the case of 
DM decays, the long lifetime results partly from the symmetry that forces one Higgs to be leptonic. We 
conclude in section IVIIi 

II. THE HIGGS AS A MESSENGER OF DARK MATTER SIGNALS 




FIG. 1: Three sectors and their interactions: the known quarks and leptons and their SU{3) x SU{2) x U{1) gauge 
interactions, a Higgs sector and a WIMP Dark Matter sector. 

What is our best guess for the structure of particle interactions at the TeV scale? In addition to 
the known physics of quarks and leptons interacting via SU{3) x SU{2) x U (1) gauge interactions, we 
expect new physics to include both a Higgs sector, responsible for SU{2) x U{1) symmetry breaking, 
and a dark matter sector. While dark matter need not be related to the TeV scale, the WIMP idea 
is intriguing: if the dark matter particle mass is of order the TeV scale, the order of magnitude of the 
observed abundance results from thermal freezeout using simple dimensional analysis. Thus we are led 
to the three sectors of Figure 1: the known sector of quarks and leptons and their gauge interactions, 
together with the unknown sectors of the Higgs and WIMP dark matter. 

What are the interactions between these sectors? Clearly the Higgs interacts directly with the quarks 
and leptons and SU{2) x U{1) gauge interactions, as shown, to give the observed masses. The interactions 
of the WIMP dark matter sector with the other two sectors are more speculative. In fact, there need not 
be any; the WIMP could annihilate to extremely light particles in its own sector. However, the WIMP 
idea is that the mass scales of the dark matter and Higgs sectors are related, and this strongly suggests 
some connection between these two sectors, as shown in Figure 1. Indeed, it could be that the Higgs 
and dark matter sectors are so closely connected that they merge. We assume that any direct couplings 
between the WIMP sector and the quark and lepton sector are subdominant, so that we are led to the 
pattern of connections of Figure 1. The Higgs sector is seen to be the messenger that makes the WIMP 
visible to us. The implications of this, for cosmic ray signatures of dark matter is clear. The WIMP, 
X, will be observed via annihilations or decays through the Higgs particles. A simple realization of this 
idea is to introduce a singlet scalar to the SM, since the only gauge invariant renormalizable operator is 
a coupling to the Higgs [19|] . This can be extended to super symmetric theories by adding gauge singlet 
dark matter to the NMSSM 0- 

What are the implications of the cosmic ray signals for this picture? The exciting thing is that the 
cosmic ray signals may shed light on the structure of the Higgs sector. The signals will depend on the 
nature and interactions of the Higgs states. The ATIC data suggests that the WIMP, will be too 
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heavy to be observed at the LHC, so that the colhder signals that could demonstrate consistency of the 
picture are those of the Higgs sector. 

Could the Higgs sector simply be the single Higgs doublet of the Standard Model? Depending on its 
mass, the Higgs would decay dominantly to pairs of top quarks, W bosons, or bottom quarks. In either 
case, for any mass of the WIMP, the signal in e"*" + e~ is smooth, and does not give the peak shown by the 



balloon experiments [2]|. Hence the Standard Model Higgs can lead to the PAMELA positron signal, but 
not the ATIC e"*" + e~ peak. If the latter is ignored, there is still the issue of whether dominant hadronic 
decays will produce a p flux larger than seen by PAMELA. For m^. < 1 TeV the p flux is apparenty one 
to two orders of magnitude too large. However, the uncertainties in the p signal are certainly an order 
of magnitude 2^ and, by going to larger values of m^, the p spectrum can be shifted to larger energies 
where there is no data. For decays to top or bottom quark pairs, > 1 TeV is in any case needed to 
explain the PAMELA signal. Values of around 100-200 GeV are possible for decays to W pairs, but 
this leads to some tension with data from anti-protons and gamma rays and also requires that the energy 



loss rate for positrons be larger [2l|, |23|. 



In this paper we concentrate on explanations of both the PAMELA e"*" data and the ATIC + e~ 
peak. Quite generally this requires dominant cascades directly to charged leptons [21'], so that we are 
immediately led to a Higgs that couples predominantly to leptons. Such Higgs bosons have received 
very little attention since they are not immediate consequences of either super symmetric or grand unified 
theories (although this is possible with some model-building). For earlier work on Higgs bosons coupling 



to leptons with different motivation, see 2j]. However, a leptonic higgs is the most straightforward 
implication of assuming that the cosmic ray data is explained by the annihilation or decay of WIMPs, 
through Higgs messengers, as illustrated in Figure 1. Furthermore, the absence, to very high accuracy, 
of fiavor violation in the charged lepton sector suggests that there is a single leptonic Higgs doublet. Hi. 
Thus, the experimental consequences for cosmic ray observations result from the connections below, while 




those for the LHC involve production of leptonic Higgs states followed by their decay to tau leptons. 

The ATIC peak in the e"*" + e~ channel will be broader than in theories where the signal results from 
X cascading directly to e or ^u. Also, there is the possibility of an energetic gamma-ray and neutrino 
component of the cosmic rays component in the energy range 100-1000 GeV. Both these signals could be 
detected in future experiments like GLAST/FERMI, VERITAS4 (gamma rays) and Hyper-Kamiokande, 
ANTARES, KM3neT, ICE-CUBE (neutrinos). These latter signals depend crucially on whether the 
cascade above results from the annihilation or decay of halo dark matter. We consider both possibilities 
in this paper. 

Consider first the case that the cosmic ray signals arise from the annihilation of dark matter in the 
halo to Hi states. The size of the signal requires a galactic annihilation cross section that is significantly 
larger than the annihilation cross section required for a successful thermal freezeout abundance, so that 
a non-thermal production mechanism 25|] is necessary. In the absence of thermal freeezeout one might 
wonder whether the motivation for WIMP dark matter is lost. Clearly the answer is no: the cosmic ray 
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signals have directly measured the annihilation cross section, and its order of magnitude is consistent 
with a weak scale mass - the WIMP motivation is actually strengthened. 

On the other hand, if the cosmic ray signals result from decays of the dark matter in the halo, then 
the abundance of the dark matter could be given by the conventional thermal freezeout of WIMPs. New 
physics is now needed to induce a dark matter lifetime of order 10^^ seconds. In the limit of stability, the 
theory possesses two discrete symmetries: one that prevents Hi from coupling to quarks and another that 
ensures the stability of x- The decay chain x ^ Hi involves the breaking of both discrete symmetries. If 
the discrete symmetries are spontaneously broken at the weak scale v, then the dimensionless symmetry 
breaking parameters may be of order v/M, giving a x decay rate of order ~ v^/M^, which is the 
desired rate for M ~ 10^^ GeV. While this is only a very rough order of magnitude estimate, it can be 
considered to be an extension of the WIMP idea, that the dark matter particle mass and interactions are 
governed by the weak scale. 

In section I VII we will introduce some explicit models for the dark matter sector and its coupling to 
the leptonic Higgs doublet. Each model can be considered in the "annihilation mode" or the "decay 
mode" depending on the absence or presence of the higher dimensional interactions that induce decay. 
Since these symmetry breaking effects are extremely small they will not affect the LHC signals of the 
model. Indeed, the main effect is on the size of the high energy photon and neutrino signals in the cosmic 
rays. For annihilations this depends on the square of the dark matter density, while for decays it is only 
linear in the density. The photon and neutrino signals have a much larger support from the center of the 
galaxy, or from satellite galaxies, while the lepton signal, because of propagation effects, has support from 
regions of the halo close to us. Thus it is likely to be the photon and neutrino signals that distinguish 
annihilations from decays. Moreover, as we will see, the neutrino signals tend to be more robust than 
the photon signals. 

III. THE LEPTONIC HIGGS BOSON 

Here we discuss the crucial features of the leptonic Higgs boson relevant for the PAMELA data; we 
postpone a discussion of the Higgs potential until section V. Consider a two Higgs doublet model with a 
symmetry that forces one Higgs doublet. Hi, to couple to the charged lepton sector and another, Hq, to 
couple to quarks 

Cyuk = y'i Qi u'jHl + y'J Qi d]Hq + y'i U e]Hi + h.c. (1) 

where Qi,u\,d\,Li,e\ are the quark and lepton fields with being the family index, and y^^^ are 
the Yukawa coupling matrices. 

We assume that the primary cosmic ray positron spectrum arises from cascade chains of the form 

^ ii{A) ^ fr(rT) ^ e+ + . . . (2) 

where x is some neutral dark matter particle, and H and A are a scalar and pseudoscalar of the Higgs 
sector, respectively. The chain could arise from either annihilation or decay of the dark matter particles, 
and the relative strength of the chains via H and A can vary. To suppress primary cosmic ray p, H and 
A must lie dominantly in the leptonic Higgs doublet Hi . Writing the neutral component of the two Higgs 



7 



doublets as 



_ ^ , hi + iai rrO _ „ , hq + iaq 

the i^, A and /i states may be written as 

H = cos ahi + sin a /i^; yl = cos (5ai — sin /3 o^; /i = sin ahi — cos a /i^ (4) 

where 

tan/3 = — < 1; + t;^ = ^2 ^ (174 GeV)^ (5) 

and the mixing angle a, that diagonalizes the Higgs mass matrix, is also taken to be small, 

sin a <C 1. (6) 

While A is the only pseudoscalar, we call the scalar orthogonal to the higgs boson H as h. Since it 
lies dominantly in Hq, and since Vq ^ vi, it is the scalar most closely related to electroweak symmetry 
breaking. It has many of the properties of the Standard Model Higgs boson, except that its couplings to 
leptons are not standard. Also, precision electroweak data do not require that rrih is close to the present 
experimental bound, since the other scalar states, and possibly the dark matter sector, may contribute 
to the S and T observables. In addition there is a charged scalar, H'^, that lies dominantly in Hi. The 
couplings of H, A, h and to quarks and leptons are proportional to the diagonal quark and lepton 
mass matrices, mu,d,e'- 

if /sin a cos a, „, 

Lh = —7^ ^{umuu +dmdd ) + -^—[emee ) 

^/2v \cosp smp 

Ca = " r- ( ■^{umuu'^ + dnidd'^) + ^{erriee'^ 

sin p 

£h = r- { —{urriuu'^ + dmdd^) H — — -(emec'^ 



iA 


( — sin /3 




cos (5 


h 


( — cos Q 




cos (3 


-( 


' — sin /? 


V \ 


cos /3 



sin /3 
sin a 
sin/? 

j-j- tjiij. I COS /3 \ 

Lu± = — ( "^^T {uVcKMmdd" + -u"^' niud^) + -^—-{uemee") (7) 
V \ COS p smp J 

where Vqkm is the CKM mixing matrix of the charged current quark interaction. 

The LEP experiments have placed bounds on the scalars H, A and . The cross section for the 
process e^e~ — > ZH is proportional to s\v?{a — j3) and depends on rriH. The limits on siv?{a — j3) are 
shown in Figure 2c of [26] for the case that H decays dominantly to rr, as we expect. For sin(a — /?) = 0.3 
the H mass cannot lie in the region (30 - 100) GeV, but for sin(a — (3) < 0.2 there is no limit. More 
importantly, the LEP experiments have searched for the process e^e~ HA which is proportional to 
cos^(a — /3) ~ 1. In the case that both H and A decay dominantly to rr as we have, the limits are shown 



in Figure 4d of 26|] and require that either mn + ttia < 20 GeV, which is strongly excluded by the width 



of the Z boson, or 



ruH + ruA > 185 GeV. 



(8) 
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If the charged Higgs, , is hghter than the th and WZ thresholds, it will dominantly decay to rvr, since 
the cs final state has a relative supression of tan^ /3. Thus the limit set by the ALEPH collaboration is 



27| 



mH+ > 88GeV. (9) 

In addition to these LEP constraints on the higgs masses, we require that the H, A ^ tt branching 
ratios are dominant and at least 0.9. This is because the dark matter annihilation or decay chain passes 
through H/A states and we need to satisfy the p constraint from PAMELA. In particular to avoid final 
states involving electroweak gauge bosons we impose 

\mH — ttiaI < mz (10) 

in order to suppress the mode H ^ AZ {A ^ H Z), and 

niH < 2mw- (11) 

to avoid H 2W. Finally, we require 

niH < 2 mA (12) 

in order to forbid the decay H — > AA since such a cascade of H would lead to a less prominant peak 
in the lepton cosmic ray spectrum for ATIC. To summarize, these constraints impose limits on mn and 
mA- 

2mz 



3 

mz 
2 



< mA < 2mw + mz 

< mn < 2m\Y (13) 



Under these conditions, H and A predominantly go to fr. The decay to hb is suppressed due to the 
leptophilic nature of H and A. 

How small should sin a and sin (3 be taken in order that the H and A couplings to quarks are small 
enough to sufficiently suppress the p cosmic ray flux? For the mass range we are considering in this 
paper, the dominant decays are to bh and rr. In this case, the ratio of quarks to leptons in the decays of 
H and A is 

i 

ml 

and 



rf = 3 ^ tan^ a tan^ (j (14) 



2 

r^ = 3^tan4/3. (15) 
^ mi; 

There is considerable uncertainty in the limit that the PAMELA p data imposes on r^. For example, a 
limit of < 0.1 can be satisfied by taking sin a and sin/3 both < 0.25. The limits on the mixing angles 
become much more stringent for the case of heavier H and A, where decays to gauge boson or tt are 
possible. In this case, sin a and sin f3 has to be sufficiently small. Although the required small mixing 
angles remain within acceptable value (the r yukawa coupling does not become large), for simplicity and 
concreteness we will not consider this case for the rest of this paper. 
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IV. ASTROPHYSICS SIGNALS 

DM particles in our galaxy and neighboring galaxies can annihilate or decay into Standard Model 
(SM) particles leading to production of cosmic rays such as electrons and positrons, protons and anti- 
protons, photons and neutrinos, which could be observed at the earth. Therefore, an observation of these 
cosmic rays consistent with DM annihilation or decay could serve as indirect detection of Dark Matter. 
However, a given framework for DM trying to explain a signal from one set of experiments must also 
respect bounds set from all other experiments. 

In this section, we study the implications of the above framework for cosmic ray signals - positrons and 
photons in particular. We also make comments about implications for the cosmic neutrino flux from the 
Galactic Center (GC) at the end. In order to carry out the analysis, we make the assumption, motivated 
in the previous sections, that DM particles dominantly annihilate or decay into higgs particles which 
have dominant Hi components. This can be guaranteed by postulating a symmetry (either discrete or 
continuous)^ and mass range for H and A given in the previous section. 

We see that the robust and distinctive feature of the above framework is production of tau leptons 
which subsequently give rise to cosmic ray electrons and positrons, as well as photons and neutrinos. 
However, the precise signal for cosmic rays depends on model-dependent details. These can be broadly 
classified into five classes. For annihilating DM, be it bosonic or fermionic, one has the following: 

XX HA {HH, AA) tttt (16) 

For bosonic annihilating DM, it is also possible to have: 

XX ^ H+H- ^ fTUrl^r (17) 

whihc is forbidden at s-wave for a majorana fermion DM by CP conservation. On the other hand, for 
decaying DM, the signals for fermionic and bosonic DM are different. For bosonic DM, there are two 
possibilities: 

X HA{HH, AA) TTTT 

X ^ TT (18) 

The second possibility in (jlSp can arise if the DM particle mixes with H or A which decays to r r. Finally, 
for fermionic decaying DM, one has the possibilities: 

X ^ H {A) V — > TTV 

X H^l^ ^ T^l^Ur] l = e,^i,T (19) 

In section we will construct some simple models which exhibit all the above possibilities. Although 
a wide variety of signals for cosmic rays can arise within this framework, it is important to note that 
the annihilation and decay modes can be related. More precisely, for the same given state (for example, 
the 4 r state in (jl6p and (jlSp ). the signal for positrons and electrons in the annihilation mode for a DM 



^ Please see section IVTl for some explicit models. 
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particle with mass and cross-section (cr v) corresponds to that for a decaying DM particle with mass 
2 and lifetime given by: 

(20) 



Ps {crv) 



where ps is the dimensionful constant appearing in the DM profiles^. The above holds true to a very good 
approximation since the electrons and positrons observed at the earth come from a short distance in the 
galaxy^ where differences in the various DM profiles are not important. Thus, the different dependence 
on the DM density profile (~ p for decays versus ~ p^ for annihilations) does not have a big effect. 



A. Positrons (& Electrons) 



In this subsection, we will estimate the positron fraction (to be compared with PAMELA) and the 
total flux of electrons and positrons (to be compared with ATIC) as a function of the mass of the DM 
and that of the higgs particles niH and mA, consistent with the assumptions above. For concreteness, 
we will show the results for the annihilation channel (I16p in which DM annihilates to 4 r's via two 
intermediate higgs particles A. This can be easily translated to results for the first decay mode in (jlSp 
from arguments mentioned above. Similarly, results for the second decay mode in (|18p can be translated 
from that obtained for the direct annihilation mode XX ~^ ^'^ i^i [^ij • We will review this result and also 
comment on the annihilation mode in (jl7p and the decay modes in (jl9p at the end of the subsection. 

The cosmic-ray background of nuclei and electrons is believed to originate from supernovae remnants 
but is not fully understood. The nuclei and electron spectra is assumed to arise from an injected flux 
which follows a power law as a function of energy, and is then propagated through the galaxy within 
some "propagation models". In the course of propagation through the galactic medium, a secondary 
component of electrons and positrons is generated by spallation of the cosmic rays on the interstellar 
medium. The parameters of the source spectra and the propagation models are constrained by fitting to 
astrophysical data. For example, the nuclei source spectra and propagation parameters are constrained 
by fitting to the proton data, the Boron-to-Carbon (B/C) ratio and so on. Since background positrons 
are dominantly generated from spallation of nuclei, this constrains the background positron flux as well. 
The electron source spectra is mostly constrained from experiments measuring the total electron flux. 
Thus, in the absence of a complete theoretical understanding of the processes involving the production 
and propagation of these cosmic rays, there is a considerable amount of uncertainty in the background 
electron and positron flux arising both from uncertainties in the nuclei and electron source spectra, 
production cross-sections, energy losses, as well as those from parameters in the various propagation 
models. It turns out that the uncertainty in the background electron spectrum is larger than that in the 
propagation parameters at present. It is important to keep these facts in mind when one tries to explain 
the data observed by PAMELA and ATIC. 

In addition to uncertainties in the background flux, there also exist uncertainties in the "signal" 
component assumed to arise from the annihilation or decay of DM particles. Once the injection spectrum 



^ For a given profile, ps is constrained by requiring that py\j = 8.5 kpc) ~ 0.3GeV/cm^. 
^ more on this in section HV Al 
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of positrons at the source is specified, the primary positron flux ^^""^ at the solar system arising from 
DM annihilation in the Milky Way galactic halo is found by solving a diffusion equation (with cylindrical 
boundary conditions for a cylinder of half-height L=l-15 kpc and radius R=20 kpc) with a source function 
given by: 



^2 , 



X 

where '^^^t {E') is the injection spectrum of positrons produced from DM annihilations or decay, and 
Pxi^) is the density profile of DM in our Galaxy. As for the background, uncertainties exist in the 
propagation parameters. Some of the most important parameters include the half-height of the diffusion 
cylinder L, the parameters characterizing the diffusion process - K = Kq f3TZ^ , where Kq is the diffusion 
constant, (5 is the velocity of the particle and TZ is its rigidity, defined as 7^ = \p{GeY)\/Z with Z as the 
atomic number, and the characteristic time for energy loss te- Different sets of parameters {L,Kq,5} 
exist which are consistent with astronomical data such as the B/C ratio, etc. The energy loss time te has 
an uncertainty of about a factor of 2 [28]. The dependence on the DM profile is weak since the positrons 
come from a short distance (O(lkpc)) where the different DM profiles are quite similar. However, the 
average local DM density, i.e. Px{r = 8.5 kpc) is itself uncertain by a factor of 2 [29^1. From (j2ip . we see 
that this uncertainty can be accommodated by a simple rescaling of the cross-section (or decay width). 

Keeping the above facts about the signal and background fluxes in mind, we have estimated the 
positron fraction and the total flux of electrons and positrons. For concreteness, we have used the 
parameters for the MED propagation model for the background and signal fluxes, as defined in 30(]. The 
background electron spectral index a defined by ^^^J^ ~ is taken to be 3.04, which is quite reasonable 



and is consistent with observations of all experiments [a] . The normalization of the background electron 
spectrum is determined by a similar procedure as described in Finally, we have used the Bessel 
approach to compute the Green's function of the diffusion equation. In particular, we have used an 
approximation to the Green's function found in 11] to solve the equation. The results of our analysis 
are shown in Figure [2j 

We see from Figure [2] that the results for the positron fraction and the total electron and positron 
flux are sensitive to the mass of the DM as well as the "Boost factor" for electrons and positrons 
Bf^^. The results are qualitatively consistent with that of 0] who have looked at a similar annihilation 
mode, although from a very different theoretical motivation and with much lighter masses of intermediate 
scalars. One finds that a DM particle with mass ~ TeV may explain the PAMELA data; however 
the explanation of both PAMELA and ATIC data requires that m-^ is about 4 TeV and Bf^^ ~ 10000. 
Also, it turns out that the above results have a mild dependence on the masses of the intermediate 
higgs particles mH,mA as long as 50 GeV ^ mH,mA ^ few 100 GeV, which is the expected range of 
masses for H and A. Note that since the results depend on astrophysical parameters such as the electron 
spectral index a and the propagation model parameters as explained above, it is possible to fit the data 
for different (but comparable) values of and boost factors Bf^^ by choosing a different combination 
of these parameters. Also, if the ATIC data is ignored, the best fit values of the parameters and 
Bfg^ will be slightly different than when both data sets are taken into account. Hence, the above results 
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FIG. 2: Results for the positron fraction (Top) as a function of energy for DM masses = 1 TeV (dotted), ~ 1.5 
TeV (dashed) and = 4 TeV (black) and mH,'mA ~ 100 GeV, in the annihilation mode ((16} with "boost factors" Bf^t 
given by 1200, 1950 and 10000 respectively. The resuh for £;3($*°_'"' + $t°t«') (Bottom) as a function of energy is only 
shown for — 4 TeV. The dot-dashed curve in both plots stands for the background. The background electron spectral 
index a is taken to be 3.04 and the MED propagation model 30] is used. Also, the reference values of the local DM density 
po and te are taken as 0.26 GeV/cm"^ and 10^'' seconds respectively. The boost factor Bt^t is explained below (I22|l . 

should only be taken as an estimate. We have not tried to optimize the fit (by a analysis) by the 
choice of astrophysical parameters. 

The boost factor Bf^^, which is a combination of various factors, deserves some explanation. More 
precisely, the boost factor is given by: 

^tot ~ ^crv ■ -B^iiimp ■ ^po ' (22) 
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Here, B^^ is the enhancement factor in the cross-section compared to the "standard" one {crv)std = 3 x 
10~^^cm^/s. -Bc/ump corresponds to the enhancement in the positron (electron) signal due to dumpiness 
in the DM halo. Strictly speaking, -Bginmp ^ function of energy [3l|. Over the relevant energy range 
of 1-1000 GeV, -B^;„^p of order few is reasonable. Bp^ is a possible enhancement due to a factor of two 
uncertainty in the local average DM density itself, as mentioned earlier. Note that a factor of two in pQ 
appears as a factor of four in Bp^ because the flux goes as Pq. Finally, the factor of two uncertainty in 
the energy loss time te mentioned above can be folded into a possible enhancement because the flux is 
directly proportional to te- Thus, Bf^f- ~ 10000 can arise in many ways. In particular, it is perfectly 
compatible with B(jy ^ 1000. In addition, the boost factor for neutrinos will in general be different than 
that for positrons (electrons) . All these facts will be important when we look at constraints from photons 
and neutrinos in the following subsections as well as the section on explicit models in section IVTl 

Moving on to the other annihilation and decay modes, the energy spectra for rs and vs in the anni- 
hilation mode (fTTll is expected to be roughly the same as in (fTHl) because m,- and my are both negligible 
compared to mH±. This suggests that it should be possible to fit the data for the same as for (fTBj) 
but with twice the boost factor. However, due to SU{2) invariance, we expect XX ~^ HH which 
gives rise to four r, to also contribute with the same cross section. Taking the sum of these modes to be 
the total cross section and normalizing it to the standard thermal cross section, we need a factor of 4/3 
relative to the Bf^^ required for only the four r case. This implies Bf^^. ~ 13300 for ~ 4 TeV. 

As explained earlier, the results for the annihilation mode (|16p can be translated to the first decay 
mode in (fT8|) . A DM mass of 4 TeV with Bf^^^ = 10000 corresponds to decaying DM with = 8 TeV 
and lifetime Ty given by ^: 



J^i^B^ciump-Bpo -Br,)^ 3.2 X 10^5 s {^[b^~~B;, . Br,) (23) 



^ 2ps{av)stdBlt 

To get the results for the second decay mode in (1180 . one needs to translate the results obtained for 



the direct annihilation mode xX ~^ tt' in 2M- From [2l[, one finds that the xX ~^ ^"^ mode gives a good 



fit to the PAMELA and ATIC data for « 2 TeV and Bf^^ ^ 3000. This corresponds, for the second 
decay mode in (fT8|) . to = 4 TeV and lifetime t-^ given by: 



^ 5.3 X 10^5 s {^Bl,^^^.Bp^ . Br,) (24) 

' — I 

Finally, we comment on the fermionic DM decay modes in (jl9p . [11] studied the first decay mode in 
(1191) a nd found that between 600 GeV and 1 TeV can explain the PAMELA data. It was pointed out 
in [1^ that this decay mode can also provide a good fit to both PAMELA and ATIC. The second decay 
mode in (jlOp is qualitatively different, since the lepton / is harder than the r. For I = e, n this provides 
a contribution to the spectrum which is steeper than that coming from the r, implying that it should be 
possible to fit the data with a lighter and a smaller boost factor. 

To summarize, therefore, even though the precise signal for astrophysics depends on model-dependent 
details, the robust characteristic of the framework is that an annihilating (decaying) DM particle with 
mass of 4 (2-8) TeV and Bf^t = 0(10000 - 13000) {t^ = 10^5-26 g) can explain the PAMELA and ATIC 



'Bp^ arises due to the fact that the flux goes as pg for annihilations but as po for decays. 
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results. We now study the consequences for cosmic gamma rays and neutrinos which provide non-trivial 
constraints on the allowed parameter space as well as give rise to potential signals for future experiments. 



B. Photons 



In general, any model of DM which produces a significant number of electrons and positrons in the 
local region of our galaxy in the energy range 10-1000 GeV (to explain the PAMELA and ATIC signals), 
is expected to dominate the production of electrons and positrons in the galactic center (GC) since the 
density of DM is expected to be much larger there. This will in turn give rise to a large yield of photons 
from inverse compton scattering (ICS) in the energy range 1-1000 GeV. In addition, there could be other 
mechanisms of photon production from DM annhilations or decays, such as final-state radiation (FSR) 
of photons from charged particle production, DM Bremsstrahlung, or from vr'^s produced from r decays. 

The total yield of photons is higher for annihilations than for decays since the flux <I> ~ for 
annihilations, while $ ~ /j-,^- for decays. This can be seen from the general expresssion for the differential 
photon flux in a solid angle region AO: 



(AO F 



1 {av) ^ ^dN^, 



47r 2m2 
1 1 



EM 



i Jannih 



An 



(25) 



47r 



dE^ 



)j Jdecay AO, 



Jannih = ^ [ dVtj p'i{r{s)) ds] Jdecay = ^ [ d9.[ p^{r{s)) ds 
JaU Jios JaU Jios 



Here, Jannih [Jdecay) Corresponds to the integrated squared (linear) DM density profile along the line- 



of-sight, and (^^)j is the photon spectrum coming from DM annihilations or decays for channel i with 



Jde 



( Jgnnih 
Psolar ^solar 



branching ratio bi. Because of the different parametric dependence on Oy, 

\Psolar solar 

especially for "steep" profiles like NFW, etc. This implies that the decay mode gives rise to a much 
weaker signal compared to that for the annihilation mode. On the other hand, the decay mode satisfies 
the existing constraints from various observations of gamma rays much more easily than the annihilation 
mode. Table IIVBI lists the values of Jannih and Jdecay (normalized such that they are dimensionless) for 
the GC for two qualitatively different profiles - the NFW and Isothermal profiles. 











Psolar r^olar 


NFW 
Isothermal 


« 15 X 10^ 
w 13 


« 28.9 
w 5.7 



TABLE I: Jannih and Jdecay for the NFW and Isothermal profiles of the GC in the Milky Way for standard choices 
of astrophysical parameters, as in 

|ii,[i3- 



Since the decay modes give rise to a much weaker signal, we will only show results for the annihilation 
mode, the one in (I16p in particular. All three sources mentioned above will contribute to the total 
photon yield for our framework in general. For photon energies E^ ^ 100 GeV, the contribution from 
ICS turns out to be the most important [3]. For E^ > 100 GeV, the contribution from vr" decay takes 
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over and dominates over the ICS and FSR contributions. However, since the quantitative predictions 
for ^ 100 GeV are subject to various uncertainties in the ICS signal from astrophysics, we do not 
attempt to analyze the above energy regime in this work. Having said that, GLAST/FERMI is expected 
to be quite sensitive in this energy range 3^]. So, a strong signal by GLAST/FERMI will provide very 
strong evidence for significant high energy electron (positron) production in the GC. 

For ^ 100 GeV, one can study the detectability of photon fluxes from DM annihilations or 
decays in Cerenkov detector based experiments such as VERITAS 4 which is expected to have a very 
high sensitivity in this energy range. The differential flux sensitivity of VERITAS 4 is expected to be 
~ 8 • 10""^- ~ 2 • 10"^ GeVcm^^g-igj.-! ^j^g energy range 100-1000 GeV [s^. The top plot of Figure 
[3] shows the predictions for the photon intensity for the annihilation mode in the 100-1000 GeV range, 
where vr*^ decays dominate the photon yield. From the figure, it can be seen that future experiments like 
VERITAS 4 have a very good potential of detecting these very high energy gamma rays in the annihilation 
mode, particularly for steeper profiles like NEW. However, steep profiles may lead to some tension with 
the flux observed by EGRET in the energy range 10-100 GeV. HESS has also made observations of 
gamma rays coming from the GC 331] and the Galactic Ridge (GR) [3j]. However, these observations 
are not ideal for DM observations because of the large contamination from gamma ray point sources 
as well as from molecular gas which are not well known. In addition, the effective Jnnnih relevant for 
these experiments is smaller because one has to subtract "off-source" contributions la]. So, we do not 
attempt to analyze constraints from these observations although one could presumably still place some 
conservative bounds. 

Another set of important constraints for E^ ^ 100 GeV comes from HESS observations the dwarf 
spheroidal galaxy (dSph) Sagittarius [35|] which is believed to have negligible foregrounds, and hence is 
suitable for DM observations. Observations of other nearby galaxies from MAGIC, CANGAROO and 
WHIPPLE give similar constraints, so we will just study the constraints from HESS. HESS reports an 
upper bound on the integrated gamma ray flux cj)^"^: _ g g ^ 10~^^ cm~^ s~^ for E^ > 250 GeV |35f ]. From 
this, the following upper bounds on the DM annihilation cross-section (decay width) can be derived: 



where = { dE^) (26) 



dK, 

where An = 2x 10"^ is the HESS sohd angle region. The quantity A^^ in (j26p . viz. the number of photons 
above ~ 250 GeV, can be estimated from the differential photon spectrum, as in [36]. At present, there is 
considerable amount of uncertainty in the DM density profile, which can vary from a cusped profile (with 
various allowed values of the cusps) to cored power-law profiles (with various values of the power-law 
exponent). We will look at two qualitatively different profiles, a large-core profile and an NFW profile, 
whose J (normalized such that they are dimensionless) values are listed in Table HVBl As for the GC, one 
finds that (— ^^^^^F — ) ^ f^IsMmih — j_ Therefore, we will only study the constraints for the annihilation 



^ P solar ^solar ) \ Psolar ^olar 

mode in (116p . 

The bottom plot of Figure [3] shows the upper bound on B^^^ as a function of rriy. for the NFW and 
Large-Core profiles of Sagittarius in the annihilation mode (jl6p . The area below the curves is consistent 
with the HESS observations of Sagittarius. A correct interpretation of these constraints requires some 
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FIG. 3: Top: The photon intensity spectrum from the Galactic Center in the annihilation mode (|16[) for DM masses 
mx = 1.5 TeV {B^^t = 1000) and = 4 TeV (B7„t = 3000) with mjf , mA = 100 GeV. The black curves stand for the NFW 
profile while the dashed curves stand for the Isothermal profile. The boost factor for photons Bl^^ is explained below (|27p . 
Bottom: Upper bound on -B^oti ^ computed in (|26|l . as a function of m-^ for the annihilation mode (|16p . The black curves 
stand for the NFW profile while the dashed curves stand for the Large-Core profile of Sagittarius. The parametrization 



in [36l | for is used. The shaded area corresponds to the uncertainty in extracting the photon spectrum from 
Monte-carlo simulations. 

explanation. Taking into account the possible enhancements in the signal from (j26p . one finds that the 
boost factor for photons is given by: 



^tot 



clump 



Bp, = Bit * ( 



B^ 

clump 
clump ^'''E 



(27) 
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Psolar "Tsolar 


NFW 
Large Core 


« 1 X 10^ 
« 1.4 X 10^ 


« 1.3 X 10"'' 
« 1.2 X 10-5 



TABLE II: Jannih and Jdecay for the NFW and Large-core profiles of the Sagittarius dwarf spheroidal galaxy for 
standard choices of astrophysical parameters, as in l^, 15|. 



where we have used (I22p . In general, the boost factors arising from dumpiness in the halo is different 
for electrons (positrons) and photons. This is because electrons and positrons observed at the earth 
come from a short distance (since they lose energy very quickly) whereas photons can come from very 
far away. This implies that B^^^^^ is close to the local dumpiness boost factor -6^°^^^ while B 



strongly depends on the direction of the gamma-ray source relative to the earth [37|]. So, 57,.. ™„ could 



7 

clump 



clump 

be smaller than from -B^/ump ^ factor around 1-10. Together with an uncertainty of a factor of two 
in Br^, this could give rise to B^^^ consistent with the bounds set by HESS. Note that the precise value 
of Bf^^ required to explain the PAMELA and ATIC data can itself be changed by optimizing over the 
astrophysical parameters consistent with the uncertainties. We have utilized the discrepancy between 
^ciump ^ciump Uncertainty in Br^ in the top plot of Figure [3l and have taken as 1000 

and 3000 for m^^ = 1.5 TeV and 4 TeV respectively^. On the other hand, if B^i^^p is not sufficiently 
smaller than -B^^^^^^, this would give rise to a tension between the parameters required for explaining 
PAMELA and ATIC (in the annihilation mode ()16p ) and the bounds set by the HESS observations of 
Sagittarius. This tension is more severe for the case of an NFW profile of Sagittarius (a little more than 
an order of magnitude) . 

To summarize, observations of gamma-rays from the Sagittarius dwarf galaxy by HESS may provide 
strong constraints for the annihilation mode (jl6p . primarily depending on the ratio of dumpiness boost 
factors Bji^j^p (in the Sagittarius direction) and -Bc/ump- the other hand, the decay modes in and 
P9P can easily satisfy these bounds. As is obvious, the situation is reversed as far as prospects for future 
signals are concerned. The annihilation modes (jl6p and (jl7p can give rise to observable signals from the 
GC, which could be measured by VERITAS 4. 



C. Neutrinos 

A DM candidate which annihilates or decays to /i's or r's will also giv e rise to a significant flux of 
neutrinos and can provide non-trivial constraints. It was noted in [l3, that the neutrino flux from 
the direction of the Galactic Center (GC) provides constraints for both annihilation and decay modes, 
with constraints for annihilation modes being much stronger (especially for steeper DM proflles). The 
neutrinos coming from the direction of the GC can be observed by detecting the muon flux induced by 
these neutrinos. A detector in the northern hemisphere (such as Super-Kamiokande) can detect upward 
going muons produced by neutrinos from the GC. One has to take into account that the three flavors 
of neutrinos oscillate into one another while traveling through the galaxy. As was pointed out in 



^ These turn out to lie in between the upper bounds set for the NFW and Large Core profiles, as seen in the bottom plot 
of Figure |3l 
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the observed neutrino-induced muon flux is almost independent of for annihilations (for fixed DM 
couplings, and small my.). This is because both the neutrino-nucleon cross-section and the muon range 
scale like energy while the DM annihilation signal is proportional to Thus, in contrast to the 

photon flux (see (j25p ). the neutrino-induced muon flux is proportional to the normalized second moment 
of the neutrino energy spectrum, i.e. oc (^)^ ^'^^ decays, a similar argument implies that the 

flux is proportional to m^. For heavier DM masses, the neutrino-nucleon cross-section grows less steeply 
as well as the energy loss term for muons starts becoming important, implying that the energetic muon 
flux is relatively suppressed [16]. It is also important to note that the neutrino flux from the direction 
of the GC is much less sensitive to the uncertainties in the DM profile, especially if one looks at the 
GC over a large-size cone centered at the GC. Then a large fraction of the total DM annihilation signal 
is contained within the observed region. Therefore, it is better to look at bounds set by Super-K for 
large-size cones (~ 10°— ~ 30°) around the GC for robust results. 



It was shown in 18|] that the neutrino- induced muon flux from direct DM annihilation to t'^t for 
2 TeV, and ^ 4500 and an NFW profile which provides a good fit to the PAMELA and 



ATIC data, is slightly above the upper bound set by Super-K [391] for a cone-size of about 10° around 
the GC. When the DM annihilates via higgs messengers to r^r", the neutrinos are softer than in the 
previous case. However, since the electrons (positrons) are also softer, explaining the ATIC data requires 
that the DM mass in this framework is quite heavy (~ 4 TeV for the annihilation mode) with a large 
)oost factor Bf^f. (« 10000). This effect tends to compensate the effect of soft neutrinos mentioned above 



la ] and one expects to get approximately the same neutrino flux as for the direct annihilation case^. 
This naively implies that the neutrino-induced muon flux for the annihilation mode (jl6p with an NFW 
profile of the GC is also above the bound set by Super-K. However, as for gamma-rays, one has to keep 
in mind that B^^^ ^ Bf^^ in general. The dumpiness boost factor -B^^^^^ is similar to that for photons 
because neutrinos, like photons, hardly lose energy and come from far away. Therefore, B'^^^^^ is strongly 
direction-dependent in general. To get more precise constraints from neutrinos coming from the direction 
of the GC, one has to know the ratio of -B^iMmp (™ direction of the GC)^ and Bf^^. As for photons, 
the results for the annihilation mode (jl6p for the isothermal profile and for all the decay modes (for most 
profiles) are within the bounds set by Super-K. 

It is worth commenting about future experiments which are going to look at the neutrino flux coming 
from various sources. The Hyper-Kamiokande experiment, which is expected to have a sensitivity bigger 
than Super-Kamiokande by about two orders of magnitude, should be able to robustly detect a positive 
signal from DM annihilations (and even decays in many cases) to neutrinos through hig gs messengers. 



For small cone-sizes (< 2°), future experiments like ANTARES 40(] and KM3neT 41[ show exciting 
prospects for the framework. For example, ANTARES and KM3neT should be able to easily observe a 
significant neutrino signal from the GC for the annihilation modes in (I16p and ()17p for steep DM profiles 
(like NFW) because Jannih increases rapidly for small cone-sizes. The prospects for the decay modes are 
not as promising since Jdecay increases very slowly for small cone-sizes. For high-energy neutrinos (^ TeV) 
arising within the framework, KM3neT may also be able to identify tau neutrinos which would greatly 
help in suppressing the atmospheric background since there are a lot fewer atmospheric tau neutrinos. 



® This has to be confirmed by exphcit analysis. 
It is expected that B^i^^p (GC direction) < -Bciump ~ B]!^!^'^^ [stI ] which would help in relaxing the bounds. 



ICE-CUBE 
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42l | does not look toward the GC, but will instead look for a neutrino- induced muon 



flux arising from DM annihilation or decays of DM particles which accrete in the earth and in the 



sun. 



43] has studied the prospects for such a flux for direct DM annihilations to neutrinos and also via 
annihilation to charged leptons such as taus. It was found that direct DM annihilation to monochromatic 
neutrinos has good prospects for ICE-CUBE. This implies that the production of neutrinos via cascade 
decays through higgs messengers, as for modes p^ . (fT7|) . ([18]) and the second mode in p^ . are also not 
promising. The first decay mode in (jlOp . however, does lead to a monochromatic neutrino, so one would 



expect that this provides much better prospects. As seen from Figure (2) in [43|], discovery is possible 
for monochromatic neutrinos if the enhancement factor in the annihilation cross-section compared to 
{o'v)std is -Bo-t; ^ 100-1000 because then the earth reaches equilibrium by the present time. However, 
the difference in the leptonic higgs framework is that the monochromatic neutrino signal arises from the 
decay of DM particles rather than their annihilation. Since the denisty of DM particles inside the earth 
p^rth -g jj^^g]^ larger than that in the galactic halo and since the flux in the decay mode scales as p^^^^ 
in contrast to as (p^"''*'^)^ for annihilations, the flux is greatly reduced [43] and hence does not provide 
good detection prospects at ICE-CUBE ^. 

To summarize, observation of neutrinos originating from DM annihilations (or decays) in the future 
will be crucial in greatly strengthening the DM interpretation of the PAMELA and ATIC signals over 
conventional astrophysical sources like pulsars since those do not give rise to a large flux of neutrinos. 
The annihilation modes and p!7|) provide stronger constraints, but also provide a greater potential 
for detectability in future experiments. Therefore, these deserve more detailed studies. It is interesting to 
note that within the decay mode, neutrinos provide a better opportunity for future detection of heavier 
DM compared to photons since for a given DM density profile and decay width, oc uiy^ for neutrinos 
which is not true for the case of cosmic-ray photons. 



V. HIGGS PHYSICS - COLLIDER SIGNALS 

We now move on to studying the Higgs sector in greater detail and potential signals for the LHC. 
In order to do that, it is important to study the higgs potential which is relevant for understanding 
the production and decay modes of the various higgs bosons. As mentioned earlier, we work within the 
framework of a CP-invariant two-higgs doublet model. After electroweak symmetry breaking, this gives 
rise to two CP-even higgs scalars h and H, a CP-odd higgs scalar A, and a charged higgs scalar . The 
couplings of these scalars to fermions was already discussed in section IIIIl In the following subsections, 
we first study the higgs potential within the framework of a leptonic higgs and then discuss signals at 
the LHC. 



* This is true even if one assumes the "best-case" scenario that the earth has reached equilibrium at the present time due 
to a sufficiently large annihilation cross-section. 
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A. The Higgs Potential 

As stated in section IIIH we consider a two higgs doublet model in which a symmetry forces one of the 
higgs doublets, Hi, to couple only to leptons, and the other higgs doublet, Hg, to couple only to quarks. 
A simple example of such a symmetry is a discrete Z2 parity. Pi, under which Hi is odd, Hq is even, the 
left-handed leptons are odd, while the left-handed quarks and right-handed quarks and leptons are even. 
Such an assignment enforces the couplings mentioned above. 

The most general CP-invariant two higgs doublet potential consistent with the above parity can be 
written as: 

V = -filiHlHq) + t^KHjHi) + ^XiiHlHqf + ^X2iHj Hif + UH\Hi){HlHq) + \,{Hj Hg){HlHi) + 

{^X5{HjH,f + h.c.). (28) 

It is not fine tuned to have a vacuum where the vev of the Hi is smaller than that of Hq by, say, a factor 
of three. A larger hierarchy of vevs can be naturally obtained as follows. The above potential has an 
asymmetric phase where Hq acquires a vev v = 174 GeV while Hi has no vev. This phase of the two 



Higgs doublet potential was studied in [45l | for the Inert Higgs Doublet model, where it was found that 
this phase has a parameter space of comparable size to the standard phase, and depends essentially on 
the sign of fif. Suppose one now introduces a small, soft, parity breaking interaction in the potential 

AV = -ifi^HlHi + h.c). (29) 

Inserting the vev of Hq into this interaction generates a linear term in Hi, and therefore a vev for Hi 
proportional to the small symmetry breaking parameter /i^, which can naturally be taken as small as 
desired. This will then guarantee that tan (3 and sin a are suppressed, as required from the arguments in 
section Hill It is convenient to parameterize the small parameter /i^ as = 2€v^ for later use. It is also 
helpful to list the number of independent parameters. Equations (|28p and (j29p have eight parameters. 
However, electroweak symmetry breaking gives rise to one condition among these parameters, reducing 
the number of independent parameters to seven. These can be taken as {Ai, A2, A3, A4, A5, e, tj^}. In the 
linear approximation for sin a and tan P (valid since both are small) , the physical higgs masses and the 
higgs mixing angle, sina, can be computed in terms of these parameters as: 

ml « 2{et-^)v^; ml ^ 2{\i) ; (30) 
m\ ^ 2{et-^ -\^)v^; m]j± ^ 2{€tf - + \^))v^; 
(A3 + A4 + A5-et;^^) 

where Sa = sina and tp = tan/3, and a and /3 lie in the range — 7r/2 < a < it/2; < /3 < ^. For 
-^1,2. .,5 = 0(1)) and /x^ = ev'^ <^ ~ ~ v'^, the parameter e = 0{l)tp. From (f30l) . this implies 
that under these conditions all higgs masses are comparable to each other, up to factors of 0(1). Also, 
depending on whether (et^^) is smaller or greater than Ai, mn could be lighter or heavier than m^. As 
we will show below, various choices of these 0(1) numbers, consistent with the constraints from LEP as 
reviewed at the end of section Hill can lead to a rich phenomenology at the LHC. 



21 



B. Potential Signals at the LHC 

The leptonic higgs has very interesting colhder phenomenology, with distinctive signal characteristics 
that hold irrespective of the explicit DM model. If the ATIC data is confirmed, then in our scheme the 
mass of the dark matter is too large for it to be made at the LHC. Rather the LHC signals are encoded 
in the Higgs messengers, and cover a wide range of dark matter models. Although Hq and Hi mix, the 
mixing is required to be small from dark matter considerations, so that the characteristics of the leptonic 
higgs is kept intact. In particular, the couplings of the mass eigenstate higgs bosons are shown in ([7]). 
Since a and (3 are both necessarily small, the higgs associated with EWSB, h, has SM couplings to quarks 
but couplings to leptons that can be considerably larger or smaller than in the SM. In the leptonic Higgs 
sector, H, A and all have enhanced couplings to leptons and suppressed couplings to quarks. This 
pattern of couplings is quite unlike that of the MSSM where the fundamental distinction is between up 
and down/lepton couplings, rather than between up/down and lepton couplings. The Higgs sector is 
as rich as in the MSSM, and hence here we are only able to provide a limited survey of the interesting 
signals. We choose to highlight the decays of the neutral higgs bosons to pairs of r leptons, as this is 
the most direct link between cosmic-ray and LHC signals. In particular, we consider each of the cases: 
Z* HA — > 4r; H, h ^ 2t and h 2A/2H 4t, in some detail and also briefly mention other 
possible LHC higgs signals such as the 8r signal and the charged higgs signal. 

1. The Z* — > HA — > tttt Signal 

The LHC signal for tttt via HA production is particularly robust because the HA production cross 
section is insensitive to the mixing angles a and /?, as they are both small, and H, A decay to rr with 
branching ratios larger than about 0.9. The cross section for HA production at the LHC is dominated 
by s-channel Z exchange (Drell-Yan production) [47'] . This is shown in Figure HI 

For example, with niH = 100 GeV and niA = 80 GeV the cross section is about 350 fb, while with 
niH = mriA = 160 GeV, the cross section is about 50 fb. This is the same as the cross section for HA 
production in the decoupling limit of the MSSM. However, in the MSSM the fr branching ratio of H 
and A is never larger than 0.08, and is frequently significantly smaller, depleting the tttt signal by at 
least two orders of magnitude, whereas for the examples of H and A masses given above we expect about 
10500 and 1500 tttt events (before imposing any cuts) with an integrated luminosity of 30 fb~^. The 
main background for this signal is expected to be diboson {Z) production in which the Z's decay to rr. 
Given that the signal is independent of the mixing angles a and /?, and depends only on mn and m^, 
which we expect to be bounded by about 2Mw and mn + mz respectively from earlier cosiderations, 
it will be very important to study a realistic simulation of this signal with selection cuts and detector 
effects. 

Before moving on, it is worth pointing out that the same Drell-Yan process will also lead to H^H~ 
production with roughly the same cross-section; hence it could give rise to an observable signal with 
sufficient luminosity for the dominant H~^ — > t^Vt and H^ t^Vt channels. This suggests a completely 
different strategy to search for charged higgs bosons in contrast to that for the MSSM, which focusses 
on production of charged higgs boson production in association with top quarks and then studying their 
hadronic (tb) and leptonic {t'^v) decay signatures. 
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FIG. 4: Cross-section for pp ^ Z* H A ^ At as a, function of ruH for different values of ruA at y/s — 14 TeV, computed 
with CalcHEP H- 



2. The h, H ^ tt Signal 

The two T signal from higgs decays has been studied in the SM. Although gluon fusion is the dominant 
production channel for the higgs, the search strategy for the hsM t't mode consists of exploiting the 
vector boson fusion (VBF) channel for hsM production, which has some distinct features that help to 
suppress the otherwise large backgrounds. Higgs production in this channel is usually accompanied by 
two jets in the forward region originating from the initial quarks from which vector bosons are emitted. 
Another feature is that no color is exchanged in the central hard process, leading to low jet activity in the 
central region. This is in contrast to most background processes. So, jet tagging in the forward region 
together with a veto of jet activity in the central region can help in achieving a high signal significance. 

We estimate the cfvbf x BR[fT) for h and and compare it with the SM case for the same higgs 
mass. Since A does not couple to WW , one only has to consider h and H. One finds: 

/ sin^ a \ 

avBF{h) X BR{h ^ tt) ^ [4^^ x BR{hsM ^ rr)] ^^^^ (31) 

[1 + (l^ - l) B^^^SM - tt)] 

avBF{H) X BR{H ^ tt) ^ [^f^ x BR{hsM ^ rr)] " 



BR{hsM ^ tt) 



For mHiinah ^ 150 GeV, the existing studies for the SM can be used to estimate the discovery potential 
in the proposed framework [isl. In these studies, the mode in which one of the r's decays leptonically 
while the other decays hadronically is analyzed in detail for an integrated luminosity of 30 For 
this mode, [cyy^p x BR{hsM —>■ tt —>■ Ij)], with I = ef/j. and j = jet, ranges between 45 and 155 fb 
for 145 > mhgj^,j > 115 GeV. After imposing various selection cuts to reduce the background, one gets 
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jetween ~ 4 to ~ 10 events for the signal, compared to about ~ 1.5 to ~ 3.5 events for the background 
48l |. This gives rise to a signal significance of about 3a — Aa for 30 fb^^ in the VBF production channel. 




FIG. 5: Signal Significance contours for tlie h ^ tt cliannel (left) and for the H 
luminosity of 30 fb~^. Results from Figure 14 in [49I ] have been used. 



TT channel (right) for an integrated 



For h and H, the signal is modified as in ()3ip . First, we identify a parameter region that gives 
large signals for both h and H. A large value of sin a relative to sin/3 is required, while both angles 
must be small to account for the DM signals. As a benchmark, one could take sin a ~ a = 0.35 and 
sin/3 /3 = 0.05 giving rise to 7, and sin(a — /?) 0.3. From (pOj) . such a value of can be 

obtained by choosing 0(1) numbers for the A's and e. For the same range of masses for the higgses h, H 
as used for hsM in [48[], one finds for the benchmark values: 

20) [a^^F X BR{h, 



<TvBF{h) X BR{h tt) 



(10 



I'SM 



avBF{H) X BR{H ^ tt) w (1.1 - 3) [a^j^p x BR{k 



SM 



.tt)] 
tt)]. 



(32) 



Thus, assuming that the results for the number of signal events scale in a simple way, one gets a huge 
signal significance of ~ 40(7 to ~ 60 o" for the h mode, and a significance of ~ 4. 4a to ~ 9a for the H 
mode, for 115 < mh,mi{ < 145 GeV in the VBF production channel for a luminosity of 30 fb~^. 

Of course, the benchmark values have been chosen to maximize the signal siginificance. To get a 
better idea of the allowed parameter space, in Figure [5] we show statistical significance contours in the 
ruh — (sin a/ sin /3) plane for h and in the uih — sm{a — j3) plane for H. As long as \a\ > /3 ^, while both 
are "small" (|a| ^ 0.3), it seems reasonable to claim that the tt Ij channel (especially for the h mode) 
provides an extremely robust signal and can be readily discovered at the LHC for higgs masses ^ 150 
GeV. The above analysis is valid for mh,mH ^ 2mw. For mh much greater than 2mw, the tt Ij 
channel does not seem promising as in this case h decays dominantly to WW and ZZ, which become 
the best modes to search for. 



f3 can always be chosen to lie in the first quadrant. 
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3. The h HH, AA tttt Signal 

A 4r signal is possible when rrih > mm{2myi, 2 m//}, with h A A or HH and subsequent decays 
into four r. The prospects for discovery through this mode are better when m/j ^ 130 GeV; otherwise 
the tt background becomes quite large and the branching ratio oi h ^ A or H goes down as well. So, we 
will first analyze the case 130 GeV ^ nifi > min{2mA, 2mj/}. The four r final state is expected to have 
smaller background than the four b final state. Again, the most favored production channel is VBF for 
the same reasons as in the previous subsection. The higgs-strahlung production channel with leptonic 
gauge boson decay can also provide a nice trigger and better handle on the background. However, the 
cross section for this channel is much lower than that for VBF. 

A full simulation study of the 4r channel with the r's decaying into 4 /x + Ai^r + 4i/^ is currently under 



study at ATLAS 50(] in the context of an NMSSM model in the VBF production channel of a SM-like 
higgs with the higgs decaying to AA followed by the decay to 4r's, just as in the proposed framework. It 
requires three leptons to be observed and triggers on one or two high pT leptons. CMS, on the other hand, 
is investigating the mode 4 r — > /i^/i^r? ^,T? ;^ containing two same sign muons and two r jets in the context 



of the same NMSSM model framework 50]. Although a full study is currently unavailable, some studies 
of benchmark points in the NMSSM have been performed 5l|. In these studies, the 4r channel is studied 
when two r's decay hadronically and the other two decay leptonically. For example, a benchmark NMSSM 
model, which seems to be quite similar to that for the proposed framework as far as the 4r channel is 
concerned, is given by: {rrih^ = 120 GeV, = 7GeV,BR{hi AiAi) = 0.99,BR{Ai tt) = 0.94} 
where the hi behaves as a SM-like higgs as far as coupling to quarks and gauge bosons are concerned. 
A preliminary study of this benchmark model indicates that a signal significance of about 20 a could be 



obtained for an integrated luminosity of 300/5 ^ 511 ]. 



In the proposed framework, one has a similar situation as the above NMSSM benchmark model for 
the 4t signal. The parameters ruh, BR{h — > AA) and BR(A tt) are very similar to their counterparts 
for the above NMSSM benchmark model for nih < 2mv/. One possible exception is the mass of A which 
could be much heavier than that for the NMSSM model while still being consistent with a significant 4r 
signal. For ruh = 120 GeV, one expects to get a similar statistical significance (~ 20 a) as in the above 
benchmark model for 300 fb^^. Since the proposed framework can have much heavier niA, the two r's 
from the decay of A are better separated implying that one could presumably get a better significance 
by utilizing this feature. If the signal and background are assumed to scale in a simple way, this would 
imply a statistical significance greater than 5 a even for a luminosity of 30 fb^^. Therefore, for ^ 130 
GeV, the 4t channel may be used to make a discovery. 

For heavier h, i.e. for m/j ^ 130 GeV, the above search strategy is not as promising because of 
the huge ti background, which begins to rise sharply at about Mrrrr = 140 GeV f5l|. Further studies 
from ATLAS and CMS are therefore needed to claim any significance for ^ 130 GeV. Finally, it 
is worth mentioning that in addition to studies at ATLAS and CMS, there is a study based on the 
proposed forward proton detector at the LHC, the so-called FP420 project [52]. This proposal utilizes 
the difi'ractive production pp — > pph and detects protons in the final state. The claim from this study is 
that the final state basically consists of events with no backgrounds, implying that the masses of h and 
A can be determined on an event-by-event basis. 

Finally, we would like to point out an interesting possibility arising from h pair-production. If nifi > 
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min{2 myi, 2 m//} as above, h pair-production could lead to a spectacular 8 r signal at the LHC. Higgs 
pair-production in the SM is dominated by gluon fusion [53]. For example, for a 120 GeV higgs (hsAi), 
the cross-section for hsM pair-production is about 35 fb. The cross-section for h pair production in 
our framework is the same as that for the SM. This implies that one could get about 1000 8t signal 
events (before imposing any cuts) for 30 fb-^. A gain, a detailed analysis of this channel would be quite 
interesting. 

VI. SOME SIMPLE DM MODELS WITH LEPTONIC HIGGS 

We now complete the picture in Figure [T] by explicitly constructing the DM sector and the couplings of 
this sector through the leptonic Higgs. To demonstrate the idea, we study three models: the L — L"^ — N 
model where the DM particle is Majorana fermion, and models where the DM particle is a scalar, either 
singlet or electroweak doublet. In each of these models, we consider both annihilation and decay modes. 



L-L^-N DM Model 



The dark sector in this model consists of a vector-like pair of lepton doublets L, and a sterile 
neutrino N . L has exactly the gauge charges of the lepton doublet in the SM. The model has two Z2 
parities - a chiral lepton parity Pi introduced in section IV Al under which Hi and all the lepton doublets 
L, and Li are odd while all other fields are even, and a dark parity Pjy under which only particles in 
the dark sector, L, L" and N , are odd and all other particles are even. Pd guarantees that the lightest 
particle among L, L" and N is stable and can be a dark matter candidate if it is neutral. 

This model is similar to the model proposed in 5j] except that it has an additional higgs which couples 
predominantly to leptons. The gauge couplings of L and are standard, and the other renormalizable 
interactions involving the dark sector are: 



rjiH^LN + r]2HfL''N + miLL'' + -tun N"^ 



(33) 



The nature of the DM candidate depends on the spectrum of the dark sector which in turn is determined 
by the mass parameters mi, rriN and the higgs vev {Hi) = vi. The charged components in L and 
form a Dirac fermion with mass mi- N and the neutral components of L and L" mix after electroweak 
symmetry breaking. The three neutral Majorana mass eigenstates (xi)X2)X3) and their corresponding 
masses (mi, m2, ms), are given by: 




/ 1 

0{vi/mN) 
\ 0{vi/mN) 



{Vl-V2)vi 
V2(m]v+m,i) 
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(34) 
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7713 ~ mN + 0{vi/mN) 

7712 ~ '^L + -^, tVi 

2{mN + rriL) 



2 



..2 



2(mjv - rriL) 



(35) 



li rriN > TTiL, the spectrum becomes: 

m^^j < m^± < 772^2 < '^xs (^^) 

and one gets doublet dark matter. The sphtting due to electroweak symmetry breaking guarantees that 
the hghtest Pd odd particle is always neutral, and the splitting between xi and X2 allows the model to 
evade the bound from direct detection by suppressing elastic scattering through coupling to a Z boson^°. 
On the other hand, ii mjsf < uil one gets: 

rriy^g < rriy^-^ < m^± < m^^ (37) 

and one gets singlet dark matter. Therefore, this model allows the DM particle to be either a singlet or 
electroweak doublet. It is also possible to have decaying dark matter in this model if the dark parity Pd 
is broken by a small amount. We now discuss both DM annihilations and decays within this model. For 
simplicity, we assume that only one of the modes is responsible for the cosmic-ray signals although in 
principle it is possible that both modes are comparable. 



1. Annihilating Dark Matter 

The phenomenological consequences in the annihilation mode for singlet and doublet DM are different. 
In addition to the differences in the cross-section, a doublet DM particle has an unsuppressed coupling to 
the Z boson in contrast to a singlet one which doesn't couple to the Z. For singlet DM, it turns out that 
the dominant annihilation channel in the s-wave is XX H Ahy a, t-channel exchange of the heavier 
neutral partner of L and L^. All other s-wave channels, such as ZZ, WW, Z-higgs, VF-higgs, and fermion 
pairs (//) are suppressed because of the tiny mixing between the doublet and singlet components. The 
H^H~ , HH and AA channels are also suppressed due to CP invariance. For doublet DM however, in 
addition to the xX ^ H A channel as before, annihilation to ZZ is also potentially relevant because of 
the unsuppressed coupling mentioned above. Since Z has a large hadronic branching ratio, the branching 
ratio of the doublet DM annihilating to ZZ has to be much smaller than that to i? A in order for it to 
be phenomenologically viable. The annihilation cross-sections for doublet DM are given by: 

(XX ^ HA). = J l + 0(^)) 



.(^), . . , ...... - 

167r7n2^(l + ^)2' '2mi 



'AT 

.4 ^2 



^° Elastic scattering through the higgs is also suppressed due to the suppressed couplings of the leptonic higgs to quarks. 
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From the argument above, this implies a lower bound on the coupling (ryl ~ ill) (for a given and 
m^' ) . We know from section IIVI that in order to fit the PAMELA and ATIC data, the DM mass is 
required to be around 4 TeV and the total boost factor for positrons (electrons) Bf^^ needs to be around 
10^. Taking into account the local fluctuation of the DM density and other uncertainties which could 
naturally give rise to a boost factor of about 10, this would imply an enhancement of about 10^ from the 
cross-section itself. For the doublet case, this naively implies that ^2 ~ ~ ^0- However, it turns out 
that there is a mild Sommerfeld enhancement in this model due to formation of a wimponium bound state 
by VF-exchange. This enhances the cross-section by a factor of about 10 if the mass splitting between 
Xi and is less than 0.1 GeV [13]. From eq. (l35p . this small mass splitting can be achieved with f ; < 5 
GeV corresponding to sin/3 3 x 10"^ (y^ ~ 0.3). Therefore, a perturbative cross-section ([38]) which 
is smaller than that naively required (without the Sommerfeld enhancement) by a factor ~ 10 is also 
allowed. This in turn implies that a much smaller r/^ — tjf is needed. For example, rji << 772 ~ 3.5 can 
explain the PAMELA and ATIC data. 

Moving on to the singlet case, the annihilation cross-section to i^yl is given by: 

a.(^) (XX ^HA),= ^^^3i^^av(^) (xx - H A% 

This is suppressed compared to that for the doublet case by typically a factor of 10"*^, which would require 
the yukawa couplings r/i , 772 to be much larger than the strong coupling limit to explain the observed data. 
Therefore, we do not discuss this case. 



2. Decaying Dark Matter 



It is interesting to consider this model in the decay mode since the annihilation mode is subject to 
more stringent constraints from cosmic-ray photons and neutrinos as explained in section llVi In order to 
have decaying DM in this model, the following breaking terms can be added to the Lagrangian (|33p 
consistent with all gauge symmetries: 

ACdecay = 5iiH\ UN + 52iHf Le^ + 5mH] el + S'^niLUL' + h.c. (38) 

(5's are required to be extremely small to allow the DM lifetime to be around 10^^~^^ seconds. The 
particular form of the above interaction and the smallness of 5's can be explained by symmetry arguments 
in several different ways. 

The mass mixing term in (j38p with coefficient 6[ dominates the decay process. Therefore, the parity 
Pi has to be suitably extended to give rise to sufficiently long decay lifetimes. One simple possibility is to 
extend Pi to the full chiral leptonic parity Pl^ x P/,, where Hi is odd under both, Li and N are odd under 
only Pl^, while e^, L and are odd under only Pi^. Then, the leading non-renormalizable operators are: 



{cuHjUN + 02^HlLe^ + c;,HjUe1) + c^rn^I^UU (39) 

where (pi, 4>r and 4>d are scalar fields that are odd under Pl^, Pi^ and Pd respectively, and are assumed 
to get vevs of order the electroweak scale, leading to: 
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5' values with the correct magnitude are naturally obtained for Muy ~ Mqut- For singlet DM x ^ ^ ^ 
the dominant decay is through the yukawa couplings in (j33p in which decays to a leptonic higgs and 
the dark doublet L, L'^ followed by the decay of L, L'^ by mixing with the SM leptons Lj through the 
mixing term in (j38p . Thus, the decay modes in this case are given by: 

or 

X +1^ ^ T^l^l^r (41) 

where I = e, r. This leads to different predictions for the electron and positron spectra for different 
I. As discussed in section IIV Al this channel will give a different fit from both 4r and 2t cases. In 
general, the spectra is expected to be a weighted average of the 4r, 2r, 2^ and 2e cases depending on the 
relative weight of channels with different /. This will change the fit to PAMELA and ATIC and also the 
predictions for photon and neutrino fluxes. We leave the detailed study of this case for future work. 

If the DM particle is a doublet, its dominant decay mode is by mixing with the SM lepton doublet 
Li. There is no constraint from any flavor experiments since the operators are suppressed by the GUT 
scale. So all SM leptons Li are equally likely to mix with the dark doublet. The possible decay channels 
in this case are x ~^ ^'^i^ 

W^P, H^T^. Since the vector boson channels are not only not suppressed by 
tan/3 but also enhanced by (m^/mvi/)^ for heavy DM (due to dominant longitudinal mode couplings at 
high energies), these channels always dominate even for very small tan/3 (corresponding to i/t ~ !)• This 
case is close to the model recently proposed in However, since this case does not fit in our original 
framework of the DM sector coupling to the visible sector through a leptonic higgs, we do not discuss 
this case in section IIVI 

Finally, in order to suppress terms like: 

C D {hIl.N + H];Lei) (42) 

which makes the DM particles decay dominantly to Iqq, we also assume that there is a quark parity Pq 
in the quark-ffg sector. If Pq is exact, this operator is completely forbidden. However, the /x^ term in 
the potential ([29]) breaks Pq. In order to simultaneously generate the /i^ term and suppress the above 
operator in a consistent way, we have to break the parity spontaneously by ~ electroweak scale vevs. 
One way to do this is by introducing a scalar (j)q that is odd under Pq and Pl^ . So the above operator 
is suppressed by and is suppressed compared to the dominant decay mode (the fourth term) in 

§8i). /z^ of the correct magnitude is generated if /z^ = e{4>r) {(f^q) , where e is a technically natural small 
coefficient. 



B. Singlet Scalar DM Model 

A particularly simple form of dark matter is a singlet complex scalar field The couplings of ^ 
to the SM is trivially constrained to the higgs sector, the extended higgs sector in our picture. With a 
discrete symmetry $ — the most general addition to the scalar potential is 

Ay = m| ^■f^ + A$t$(i7/i/; + xi?|//g) + A'($"^$)^ (43) 
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We assume the relative strength x of the quartic couplings to Hq and Hi is small, i.e. x < 1/3, so that 
dark matter annihilates dominantly via hJHi to At. Taking m$ = 4 TeV, to allow for an explanation 
of the ATIC as well the PAMELA data, we find a galactic annihilation cross section relative to {(Tv)std = 
3 X 10~^^cm'^s~^ by a factor Ba-v ~ 50 n^R, where A = nvr^ and R is the Sommerfeld enhancement factor. 
An attractive yukawa potential is generated between the annihilating $ particles by the exchange of the 
scalars in both Hi and Hq, with a strength proportional to X'^vf and x^A^f^, respectively. Since the 
dark matter particle is expected to be about 30 times heavier than the exchanged scalar, the resulting 
Sommerfeld enhancement factor can be significant. For example, n = 4, and either x ~ 1/3 or vi/vq ~ 1/3, 
leads to a region where R is rapidly varying from 10 to in excess of 100. Hence B^-v of order 10^ — 10^, as 
required to explain PAMELA and ATIC, is possible, even for perturbative values of A. As A is increased, 
the one-loop radiative correction to the mass term hJ Hi also increases, leading to a little hierarchy 
problem of why the H and A states are significantly lighter than the dark matter This is a general 
naturalness problem for models with annihilation of heavy dark matter to leptonic Higgs states lighter 
than 2Mw. The large coupling needed for a large annihilation cross section leads to a large radiative 
contribution to the leptonic Higgs mass parameter from a loop with internal dark matter particles. 

We can also have a decaying scalar dark matter by introducing a Zq parity so that the leading 
interaction that couples $ linearly to SM particles is the dimension 6 interaction 

AC = -^(Pl<i>HjHi + h.c. (44) 

Under the Zq symmetry, <I) — > — <i) and (pD e^^ cpm while all SM particles transform trivially. Taking 
the scalar 0d to acquire a weak scale vev, and taking muv oi order 10^^ GeV, leads to decays of <I> to 4r 
with a lifetime of order 10^^ seconds. This case would require to be around 8 TeV (see section WV\i . 



C. Inert Higgs-Doublet DM Model 



Finally, we consider the scalar DM to be an electroweak doublet. A simple example of this is the inert 



Higgs Hj, first proposed in [55|. Extending this model, we have three Higgs-doublets, Hq, Hi and Hj, 
where both Hq and Hi get a vev but the inert Higgs Hj does not. The Higgs doublets are all identical in 
the sense of gauge charges, but have different masses and parities. The couplings of Hj to fermions are 
forbidden by dark parity. Since no symmetry can forbid the term {H\Hi){HlHq), we have no symmetry 
explanation of suppressing the couplings of Hi to the quark sector in this model. One can only claim that 
it is not unreasonable to have such couplings being numerically suppressed. However, such a possibility 
is at least naturally allowed within the framework of a leptonic higgs sector. 

From now on, we concentrate on couplings to leptonic higgs. The relevant new terms in the scalar 
potential are: 

AV = m]H]Hi + H]HiH}Hi-{H]Hi){H}Hi) + {H]Hif + C.C. (45) 

We omit the 0{\) coefficients in front of each quartic term to avoid use of new notation. These terms allow 
the lightest Po odd particle to be neutral after electroweak symmetry breaking splits the doublet. The 
second quartic term splits the charged and neutral components of Hi, while the last term in (j45[) splits 
the scalar and the pseudoscalar parts of the neutral component in Hi, naturally evading bounds from 
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direct detection. The DM particles can annihilate to HJHi At. However, similar to the doublet DM in 
the LU-N model, the ZZ annihilation channel is also present. In order to fit the ATIC data we need some 
of these quartic couplings to be large and the mass splitting between the scalar and pseudoscalar (and 
charged and neutral) components to be small (to provide a modest Sommerfeld enhancement), requiring 

a hierarchy between these quartic couplings. 

A decaying scalar dark matter is also possible by introducing small dark parity breaking term. Similar 
to the LL'^N model, the lowest dimensional operator is the mass mixing term: 

A-Cdecaj/ = 5m^H\Hi + c.c. (46) 

The corresponding term involving the quark Higgs can be forbidden by imposing an exact quark parity 
as we did for the the LL^N model. The desired small value for (5m^ can again be obtained by extending 
Pi to the full chiral lepton parity and the dark parity Pd to Z^ in which Hj has 2 units of charge. Then, 
the parity preserving non-renormalizable operator is: 

^jC-decay = H\Hi + C.C. (47) 

The parities of these 0's are as defined in the LL'^N model except the <j)D — iff^D under the Z^ discrete 
dark symmetry. The possible decay channels are WW, ZZ, (A, H) h, tt and the three-body decay to 
{A,H)hh. The vector boson channels are enhanced by (m^/mw)'^ as before but unlike the case for 
fermionic DM, these vector boson channels are also suppressed by sin^ (3. Therefore, the fr channel can 
dominate when sin/? is relatively small. For example, for sin/3 ~ 10~^ (y-r ~ 1), the decay branching 
ratio to vector bosons is about 5% for 4 TeV dark matter. The scalar decay channels come from quartic 
interactions (for example A3|i7;|^|i7gp) in the potential which can not be forbidden by any symmetry. 
However, the two-body decay mode is suppressed by at least {vq/M^)'^ ~ 10~^ while the three-body 
decay mode is suppressed by an extra phase factor of ~ l/27r^. For quartic couplings not much bigger 
than unity, the rr channel can dominate the decay. Since the final state is rr, a 4 TeV decaying dark 
matter of this model is expected to fit the data well. 



VII. CONCLUSIONS AND SUMMARY 



In this paper we have explored the consequences of assuming that a leptonic Higgs mediates the 
interactions between the dark matter sector and the Standard Model sector. The motivation for this 
approach is two-fold. One is theoretical: the same TeV mass scale underlies both the breaking of weak 
interactions and the dark matter annihilation rate, indicating a close connection between the dark matter 
and Higgs sectors. The observational motivation is that the cosmic ray signals are leptonic; in particular, 
the limits on a primary p flux indicate that mediation via a Higgs with quark couplings should be sub- 
dominant. We study the minimal case of a single leptonic Higgs, in which case its largest couplings to 
matter must be to r leptons, with couplings to e and fi that are sufficiently small to play no role. The 
Higgs potential could lead to mass mixing between the states with leptonic and quark couplings but, 
again, the absence of an exotic primary cosmic ray p flux limits this mixing to be small. Hence, there is 
one neutral mass eigenstate scalar, H, one pseudoscalar. A, and one charged scalar i7+ that maintain 
their dominant leptonic couplings. This allows an important connection between the leptonic cosmic ray 
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signals and the expected Higgs signatures at the LHC. In both cases there are signals that arise from the 
production of H, A and and their subsequent decays to r leptons. From LEP limits and the leptonic 
cosmic ray signals we argue that the states H, A and are likely to have masses in the range of roughly 
100 GeV to 200 GeV. 

The leptonic cosmic ray signals, for both PAMELA and ATIC, could arise from a variety of channels: 
dark matter annihilation to r^,r^i/^ or dark matter decay to t'^,t'^ ,t^v,tvI via intermediate H,A and 
states. For concreteness we have computed the cosmic ray signals for the case of annihilations 
XX ~^ 'T^j Slid the results are shown in Figure [2j Good fits to the PAMELA and ATIC data can be 
simultaneously obtained for in the region around 4 TeV. We have also commented on all the other 
possible annihilation and decay modes. The width of the peak in the leptonic signal around 600 GeV is 
larger than in the case of annihilations or decays directly to e and ^u, and we expect this to be a common 
feature of all modes involving rs. Similarly we expect the signal in the positron fraction to continue to 
energies larger than the present PAMELA data, and not to show a sharp peak. Such an annihilation 
signal requires a large total boost factor in the annihilation cross section of order 10^, but essentially 
identical cosmic ray signals can result from dark matter decays with lifetimes of order 10^^ seconds. 

We have computed the photon spectrum that results from annihilations of DM via pairs of leptonic 
Higgs states to r'^ in the dwarf galaxy Sagittarius. For a DM mass of 4 TeV, HESS data places a limit 
of ~ 10'^ (~ 10^) on the relevant boost factor, B^^^, for the case that Sagittarius has an NEW (Large 
Core) profile. For this annihilation mode and DM mass, the PAMELA and ATIC signals require a boost 
factor Bf^^ that is of order 10^. However, these two boost factors are not identical in general, so further 
observations of hard gamma rays from Sagittarius would provide a powerful probe of this annihilation 
channel. We also compute the high energy gamma ray signal expected from DM annihilations in the 
center of the Milky Way galaxy, as shown in the top plot of Figure El which could be detected by 
VERITAS 4. The annihilation channel to 4r also produces a galactic flux of neutrinos which is in mild 
conflict with the bounds set by SuperK on upward going neutrinos through the Earth, if the relevant 
neutrino and electron boost factors are identical. As before, these boost factors can differ, so future 
neutrino measurements also offer the possibility of a crucial independent verification of the DM origin 
of the charged lepton cosmic ray signals. In addition, an observation of energetic cosmic neutrino flux 
strongly favors the dark matter interpretation of the cosmic-ray signals over astrophysical ones (such as 
pulsars) in general since astrophysical sources do not emit a large flux of high energy neutrinos. 

For the decay modes, constraints from current gamma-ray and neutrino observations are easily satisfied 
but at the same time these also lead to less promising prospects for future experiments. For heavy decaying 
dark matter, neutrinos provide a better opportunity for future experiments compared to gamma rays, 
since the neutrino flux is proportional to the dark matter mass. 

A common feaure of the cosmic ray signals discussed above is that the dark matter particle has a 
mass in excess of 1 TeV and is not expected to be produced at the LHC. While cosmic ray photon and 
neutrino fluxes allow discrimination between DM annihilation and decay, the LHC probes a completely 
complementary feature of the theory - the production and decay of the leptonic Higgs states A and 

. There is a very rich and distinctive Higgs phenomenology; all the couplings of these states to matter 
and to electroweak gauge bosons are determined by just their masses and two mixing angles, as shown 
in ([T|). The signal for HA production via a virtual Z is particularly robust because it is insensitive to 
the values of the small Higgs mixing angle a and to the ratio of vevs tan (5. Depending on the H and A 
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masses this HA production cross section is typically in the (50-250) fb region, and from the cosmic ray 
signals we know that the branching ratio to is close to unity. Detailed simulation studies are needed 
to estimate the observability of this signal. 

We have computed cjxBR for the production of both the leptonic Higgs scalar H and the Higgs boson 
h in the vector-boson fusion at the LHC followed by decay to r pairs. For mH,mh < 2M]y, 

discovery at the LHC is possible for a wide range of a and /?, provided \a\ > (3, via the channel fr Ij, 
and for some regions of parameter space the signal is very large allowing prompt discovery at the LHC. 
If rrifi > 2m H or 2mA the signal for h production is lost, but the cascade h — > HH, A A 
leads to a 4r signal instead. Comparisons with benchmark studies for 4r signals in the next-to-minimal 
supersymmetric theory suggest that a 5a discovery may be possible with an integrated luminosity of 30 
fb~^ for mil < 130 GeV. Further simulation studies are needed, especially for mh > 130 GeV where 
backgrounds from top quark pair production become important. Other LHC Higgs signals, such as 
charged higgs pair-production by Drell-Yan followed by their decay to and a spectacular 8r signal 
from h pair-production, are also possible and would be quite interesting to study further. 

Simple particle physics models of DM that couples to the visible sector via a leptonic Higgs are very 
easy to write down. In the case that DM is a heavy lepton, the interaction with the leptonic Higgs is 
via a yukawa coupling, allowing s-wave annihilations to HA. On the other hand, scalar dark matter 
can annihilate via a quartic scalar interaction to H or A pairs. In either case a large annihilation cross 
section is possible, aided in many cases by a modest Sommerfeld boost factor; but in both cases the 
DM abundance must be produced non-thermally. Alternatively, the DM could be produced thermally, 
with the cosmic ray signals arising from decays. The long lifetime follows from a decay amplitude that 
is suppressed by two powers of the ratio of the weak scale to the unified scale, which arises from a 
combination of breaking the discrete symmetry that leads to near stability of the DM and the discrete 
symmetry that ensures the leptonic nature of the Higgs. 
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